Correlated metallic state in honeycomb lattice: Orthogonal Dirac semimetal 



(N 

O 
(N 

Oh: 

00 



u 

a 

o 
o 



in 

o 

in 
o 



Yin Zhong/'Q Ke Liu,^ Yong-Qiang Wang,^ and Hong-Gang 

^ Center for Interdisciplinary Studies & Key Laboratory for Magnetism and 
Magnetic Materials of the MoE, Lanzhou University, Lanzhou 730000, Chtna 
'^Institute of Theoretical Physics, Lanzhou University, Lanzhou 730000, China 
'^Beijing Computational Science Research Center, Beijing 100084, China 
(Dated: March 4, 2013) 

A novel gapped metallic state coined orthogonal Dirac semimetal is proposed in the honeycomb 
lattice in terms ol Z2 slave-spin representation ol Hubbard model. This state corresponds to the 
disordered phase of slave-spin and has the same thermaldynamical and transport properties as usual 
Dirac semimetal but its singe-particle excitation is gapped and has nontrivial topological order due 
to the Z2 gauge structure. The quantum phase transition from this orthogonal Dirac semimetal 
to usual Dirac semimetal is described by a mean-field decoupling with complementary fluctuation 
analysis and its criticality falls into the universality class of 2-)- ID Ising model while a large anoma- 
lous dimension for the physical electron is found at quantum critical point (QCP), which could be 
considered as a fingerprint of our fractionalized theory when compared to other non-fractionalized 
approaches. As byproducts, a path integral formalism for the Z2 slave-spin representation of Hub- 
bard model is constructed and possible relations to other approaches and the sublattice pairing 
states, which has been argued to be a promising candidate for gapped spin liquid state found in 
the numerical simulation, are briefiy discussed. Additionally, when spin-orbit coupling is consid- 
ered, the instability of orthogonal Dirac semimetal to the fractionalized quantum spin Hall insulator 
(fractionalized topological insulator) is also expected. We hope the present work may be helpful for 
future studies in Z2 slave-spin theory and related non-Fermi liquid phases in honeycomb lattice. 



I. INTRODUCTION 

It is still a challenge to understand exotic quantum 
phases and its corresponding quantum criticality beyond 
the conventional Ginzburg-Landau- Wilson paradigm. i"— 
One of ideasii^"— to attack this challenging problem is 
to fractionalize the electrons in the model Hamiltonian 
into more elementary collective excitations, e.g., quasi- 
particles like spinon, holon, and so on, in terms of slave- 
particle approaches near the putative quantum critical 
points due to wild quantum fluctuations.— i^^"— 

Recently, a remarkable determinantal quantum Monte 
Carlo simulation has been performed in the honeycomb 
lattice at half-filling,— This study claimed that there ex- 
ists a short-range resonating valence-bond liquid with a 
spin gap between the Dirac semimetal and antiferromag- 
netic insulating phasesj^"— (A similar gapped spin liq- 
uid phase has also been found in the Kagome Heisenberg 
model by density-matrix-renormalization-group>22, How- 
ever, it is noted that a very recent quantum Monte Carlo 
simulation, which extend the previous study to much 
larger clusters, finds no evidence of a spin liquid region 
found in Ref. |2^] j2S) Immediately, the nature of this pos- 
sible spin liquid state has been explored by using slave- 
particle framework with U(l) and/or SU(2) gauge struc- 
ture and was proposed to be a Z2-like spin liquid^Si^i^^ 
which was carefully inspected by later numerical simu- 
lations in J1-J2 Heisenberg model. (Note, however, 
simple reduction of the physics of the Hubbard model to 
that of an effective spin-1/2 system, such as the J1-J2, 
might not be achievable in light of Ref. [s^l-)^ 

Instead of studying the elusive insulating spin liquid 
phase found in the numerical simulation, it is also in- 



teresting to see whether a possible fractionalized metallic 
state can exist in the honeycomb lattice. This is stim- 
ulated by a recent paper of Nandkishore, Metlitski and 
Senthil^ where they reinspected the Z2 slave-spin rep- 
resentation of single-band Hubbard model and pointed 
out that the correct disordered state of slave spins is not 
a putative Mott insulator but an exotic fractionalized 
metallic state called orthogonal metali^"— This state is a 
compressible metal, which has the same thermodynamics 
and transport as the usual Landau Fermi liquid, but its 
electronic spectral function has a gap, which results in a 
simplest non-Fermi liquid. Therefore, there is an orthog- 
onal metal-Fermi liquid transition instead of Mott tran- 
sition in the slave spin representation of the single-band 
Hubbard model. Furthermore, for the multi-orbital mod- 
els like the Anderson lattice model, an orbital-selective 
orthogonal metal transition has been found by the Z2 
slave-spin approach^^ and is argued to be an alternative 
Kondo breakdown mechanism for certain heavy fermion 
compounds^iii^i^sda-is 

In the present paper we try to uncover a possible frac- 
tionalized metallic state in terms of Z2 slave-spin repre- 
sentation of the Hubbard model in the honeycomb lat- 
tice at half-filling following the same methology of Nand- 
kishore, Metlitski and Senthili^ It is found that while 
the ordered state of slave-spin is identified as the usual 
Dirac semimetal, an exotic gapped semimetal named or- 
thogonal Dirac semimetal survives when the slave spin 
becomes disordered. The orthogonal Dirac semimetal 
has the same thermodynamics and transport as the usual 
Dirac semimetal but with the nontrivial topological order 
due to the local Z2 gauge structure. 

Furthermore, we analyze in detail the correspond- 
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ing quantum phase transition (QPT) between these two 
metalhc states at mean-field level and find that it falls 
into the universality class of 2+lD Ising model while 
a large anomalous dimension for the physical electron 
is found at quantum critical point(QCP). This result 
could be considered as a fingerprint of our fractionalized 
treatment when compared to other non-fractionalized ap- 
proaches. Additionally, the mean-field analysis is sta- 
ble to weak fluctuations since vanishing density of states 
of the Dirac semimetal suppresses the effect of Landau 
damping even under the condition without long-range 
Coulomb interaction. This is in contrast to the origi- 
nal case of orthogonal metal where the mean-field treat- 
ment of quantum criticality is only valid in the presence 
of the long-range Coulomb interaction. As byproducts, 
a path integral formalism for the Z2 slave-spin repre- 
sentation of Hubbard model is constructed and possi- 
ble relations to the exotic spin liquid named sublattice 
pairing states, fractionalized quantum spin Hall insula- 
tor, effective Gross-Neveu theory, the dual approach of 
interacting fermions and the gauge/gravity duality are 
briefly discussed. Moreover, we expect the orthogonal 
Dirac semimetal may be realized in future experiments 
of ultracold atoms in the honeycomb optical lattices'*^ 
and surface states of three dimensional topological in- 
sulators since they naturally support Dirac fermions as 
their low-lying excitations. ^^=52 In addition, it may be 
interesting to generalize the present analysis of orthogo- 
nal Dirac semimetal to the bilayer graphene-a semimetal 
with a non- vanishing density of state, which is known to 
have a strongly correlated ground state i^i^ 

The remainder of this paper is organized as follows. In 
Sec. ini we introduce the Hubbard model in the honey- 
comb lattice and reformulated it in terms of the Zi slave- 
spin representation. Moreover, a path integral formalism 
for this slave-spin representation of the Hubbard model 
is also constructed in this section. Then, a mean-field de- 
coupling is used and two resulting mean-field states are 
analyzed in Sec. IIIII One state is the expected orthog- 
onal Dirac semimetal and the other is the usual Dirac 
semimetal. In Sec. IIVI the QPT between these two states 
are discussed with emphasizing on the scaling behaviors 
and the stability of mean-field treatment of quantum crit- 
icality is studied. Some relations to other approaches and 
instability of the orthogonal Dirac semimetal are briefly 
discussed In Sec. [V] Finally, Sec. |Vl] is devoted to a 
conclusion which summarizes our main findings. 



II. Z2 SLAVE SPIN REPRESENTATION AND 
THE HUBBARD MODEL IN THE HONEYCOMB 
LATTICE 

The model we used is the Hubbard model in the hon- 
eycomb lattice at half-filling^i^S 

^ = E + h.c.) + ^ - 1)'' (1) 



where Ui ~ X^o- '^L'^io'' ^ onsite Coulomb energy 

between electrons on the same site and t is the hopping 
energy between nearest-neighbor sites. Since we are in- 
terested in the case of half-filling, the chemical potential 
has been set to zero. This model has been studied by 
many authors and is believed to exhibit several distinct 
phases depending on the ratio of IJ jt. (For a review, see 
Ref. For small C//i, the usual Dirac semimetal ap- 

pears with nearly free relativistic Dirac fermions being 
the low energy excitation. In contrast, an antiferromag- 
netic Mott insulator survives for large U jt and spin ro- 
tation symmetry breaks spontaneously.^^ Besides, some 
exotic spin liquid Mott insulating states with a charge 
gap, e.g. algebraic spin liquid (ASL) and Z2 spin liq- 
uids, may exist in the intermediate coupling as argued 
by slave-particle techniques with the projective symme- 
try group (PSG) analysis and inspected by sophisticated 
numerical simulationsi^i^"— 

Here, the problem we are interested in is whether 
there exists an alternative metallic state beside the triv- 
ial Dirac semimetal. It seems that such a metallic 
state is likely due to the fractionalizing scheme proposed 
by Nandkishore, Metlitski and Senthil,^^ in which they 
found a Z2 fractionalized metal named as the orthogonal 
metal in the Z2 slave-spin representation of the Hubbard 
model. The orthogonal metal was incorrectly identified 
as a non-magnetic Mott insulator by previous studies be- 
cause the slave-spin was considered to carry charge of the 
physical electron. As a matter of fact, as also pointed 
out by Nandkishore, Metlitski and Senthilj^ the slave- 
spin does not carry any quantum numbers of the physi- 
cal electron but only has the Z2 gauge charge since it is 
bound by the gapped Z2 gauge field (a brief explanation 
of this point is given by Appendix A) Thus, all quantum 
numbers of the physical electron have to be taken by 
the slave-fermion, which also carries a Z2 gauge charge 
to ensure the gauge invariance of the physical electron 
composed by a slave-spin and a slave-fermion. 

Therefore, it is expected that a fractionalized metallic 
state similar to the orthogonal metal may exist in the 
honeycomb lattice when using Z2 slave-spin representa- 
tion. To this point, a careful reader may wonder what 
the difference is between our study and the case of Nand- 
kishore, Metlitski and Senthil^^ since the same starting 
points have been made, namely, the Hubbard model and 
the Z2 slave-spin representation. We should emphasize 
that in the paper of Nandkishore, Metlitski and Senthil^ 
a square lattice is kept in mind implicitly and the refer- 
ence state is the usual Landau Fermi liquid with a large 
Fermi surface since a Fermi liquid is indeed a good start- 
ing point for a small V jt in the square lattice. How- 
ever, in the present work, we will discuss the Hubbard in 
the honeycomb lattice at half-filling and instead of the 
Fermi liquid, a Dirac semimetal is our reference state, 
which has relativistic Dirac fermion excitation near dis- 
connected Dirac points. Thus, it will lead to a different 
fractionalized metallic state from the orthogonal metal. 
This will be seen in the following discussion. Now, we 
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turn to the Z2 slave-spin framework for a honeycomb 
lattice at half-filling. 

A. Zi slave-spin representation of the Hubbard 

model 

In the treatment of Z2 slave-spin approach, the physi- 
cal electron is fractionalized into a new slave fermion 
/o- and a slave spin as^^i^ 

^ia — jia^i (2) 

with a constraint rf + I — 2{ni — 1)'^ enforced in every 
site. Under this representation, the original Hamiltonian 
can be rewritten as 

{ij)tT i 

where rn = n{ = Y^afLff^- Obviously, a Z2 local 
gauge symmetry is left in this representation (both slave- 
fermions and slave spins carrying the Z2 gauge charge) 
and the corresponding low-energy effective theory should 
respect this. The mentioned gauge structure can be 
seen if f^^ ei/fp and rf urf with u = ±1 
while the whole Hamiltonian H is invariant under this 
Z2 gauge transformation. Moreover, as argued by Nand- 
kishore, Metlitski and Senthil;^ the slave-fermion car- 
ries both charge and spin of the physical electron Cg. while 
slave spin encodes the remanent quantum coherence 
of physical electrons, which has also been mentioned on 
the issue of fractionalized quantum spin Hall insulator;^ 

B. Path integral formalism for the Z2 slave-spin 

representation of the Hubbard model 

Before turning into discuss the possible alternative 
metallic state, we present the construction of a path in- 
tegral formalism for the Z2 slave-spin approach of the 
Hubbard model in this subsection. To this aim, we fol- 
low the approach of Ref. [s^l where the general Z2 gauge 
theory is constructed in an extended Hubbard model. 

The construction of the path integral is to calculate the 
partition function Z = Tr(e^'^-^P) where P is the pro- 
jective operator to exclude unphysical states introduced 
by Z2 slave-spin representation. Here we use 

p = Y[\{i + (-i)^[-r+i-2("{-i)=i). (4) 

i 

This choice has the advantage to meet the mean-field 
theory of Z2 slave-spin approach. Obviously, one can 
employ another equivalent projective operator— 



We will use the first definition of P in the following dis- 
cussion. Follow Ref. [57|, the projective operator can be 
reformulated by introducing auxiliary Ising field cr^ = ±1 

i ^ (Ti=±l 

Since [P, H] = 0, one can define an effective Hamiltonian 
Heff as 

Heff = + E - '^Okf + 1 - 2{n{n (7) 

i 

Then using the same method in the treatment of quan- 
tum Ising model (see Appendix B) and standard coherent 
state representation of fermions, one obtains the follow- 
ing path integral formalism of Z2 slave-spin representa- 
tion of Hubbard model^ 

^ = n/ dMn)dMn)dip,{n)Si^'f-l)da,in)S{a^~l)e-^ 

ni 

(8) 

and 

S = ^fia{n)(.ha{n) - fia{n- 1)) 
-\-^^(~Pi{n)ipi{n-\- l)a(n) 

+eJ2^lil-<^^inm-2{n{)% (9) 

where we have used T[\ip) ~ ipi\ip) with ip — ±1 and 
Tflv) = \V'i)\V2)\^3) ■ ■ ■ \ - 'Pi) ■ ■ ■ I'Pn) to avoid confu- 
sion with auxiliary Ising field (7^ with a{n) — — ^ In e(^ + 
i^(l — fTj;(n))). The above action is our main result in 
this subsection and further approximations have to be 
made in order to gain some physical insights. A popu- 
lar approximation is to decouple the interaction term be- 
tween slave-spin if and slave-fermion fa at the mean field 
level and then reintroduce phase fluctuations (in fact, a 
Z2 gauge field due to the gauge structure of H)^^^ De- 
tails of this treatment has not been reported until now, 
however, we hope our formalism constructed here may 
be useful in this direction. 



III. MEAN-FIELD THEORY AND 
ORTHOGONAL DIRAC SEMIMETAL 

Undoubtedly, it is a formidable task to treat exactly 
the Hamiltonian of the Hubbard model in the slave-spin 
representation given by Eq. ^ and its path integral 
formalism given by Eq. ([8]). Thus we only consider a 
mean-field treatment here and discuss the stability of the 
mean-field analysis in the next section. 
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It is straightforward to derive a mean-field Hamilto- 
nian as foUows^S 

^/ = - E fe/l/^- + f^-^-) -'^Y.^M - i)',(io) 

{ij)<y i 

Hi = -Y. (J^^^^^! + ^■^■) + E (^^ + 7) (11) 

where the Lagrange multipher Ai has been introduced 
to fulfill the constraint on average, = t(rf rj), = 

^12a(ficrfj'r) ^^'l extra self-consistent equation ap- 
pears as {tD + I — 2((ni — 1)^) due to the constraint. The 
decoupled Hamiltonian Hj is an extended quantum Ising 
model in transverse field and Hf describes / fermions 
in the honeycomb lattice. Here, at the mean-field level, 
a further simplification can be made by setting all the 
Lagrange multiplier Xi to be zero, provided only non- 
magnetic solutions are involved and a half-filling case is 
considered.-^ This means the constraint is not violated 
seriously when magnetic order is absent. Therefore, we 
will drop the constraint term in the mean-field Hamilto- 
nian Eqs. (Uni) and pT|) hereafter. 



A. Quantum Ising model in the honeycomb lattice 

Let us first focus on the quantum Ising model given by 
Eg. (1111) . It is well known that the standard transverse 
field Ising model in one spatial dimension can be exactly 
solved by Jordan- Wigner transformation and it has two 
phases with the critical exponents being the same as two- 
dimensional classical Ising modelji Beyond one spatial 
dimension, to our knowledge, no exact solutions exist for 
the quantum Ising model until now. However, gener- 
ically, one may define (r^) as a useful order parameter 
and there are at least two phases in two space dimensions 
or beyond. (It is just this case in the study of single-band 
Hubbard model in terms of some mean-field approxima- 
tions and the Schwinger bosons theoryi^i^i^) One is a 
magnetic state with (t^) 7^ while the other is described 
by a vanished (r^) and is a disordered state with an exci- 
tation gap. Moreover, there exists a QCP between these 
two distinct phases, whose critical properties could be 
described by a quantum theory. 

In the case of the extended quantum Ising model in the 
honeycomb lattice, Riiegg and Fiete^^ found that there 
exist a ferromagnetic phase (r^) 7^ and a paramagnetic 
phase (t^) =0 by using a 4-site cluster- mean- field ap- 
proximation for the mean-field Hamiltonian. However, 
the paramagnetic phase is interpreted as the valence- 
bond solid (VBS), which turns out to break both lattice 
rotation and translation symmetry. Thus, if so, the pu- 
tative quantum critical point between these two states 
will be unstable since the ferromagnetic phase and the 
valence-bond solid breaks entirely different symmetries 
and a first-order transition is expected generally accord- 
ing to the Ginzburg-Landau- Wilson paradigm. (A de- 



confined criticalityi^ cannot be excluded in principle but 
we have not seen any signals for this consideration.) 

However, we note a different slave-spin treatment, 
which combines with a Schwinger boson analysis instead 
of a cluster-mean-field approximation, favors paramag- 
netic phases without broken lattice rotation and transla- 
tion symmetry if the onsite U is not sufficiently largei^ 
Furthermore, a slave-rotor approach has also been used 
to study the strong coupling behaviors of the Hubbard 
model in the honeycomb latticeJ^i^'^ In these papers, 
the insulating state of quantum XY model, which cor- 
responds to the paramagnetic phase of slave-spin repre- 
sentation, is found to preserve all physical symmetries. 
Thus, we expect a paramagnetic phase without any bro- 
ken symmetries in the effective quantum Ising model in 
the honeycomb lattice with the belief that the slave rotor 
and slave-spin approach are equivalent though different 
interpretations are involved Jii^i^^iiS 

Therefore, we may simply assume that the quantum 
Ising model in the honeycomb lattice has a ferromagnetic 
ordered phase ((r^) ^ 0) and a disordered paramagnetic 
phase without any broken symmetries ((r^) = 0) with 
a QCP between them. Moreover, an effective theory for 
the quantum Ising model in the honeycomb lattice can 
be derived by using its path integral formalism (see Ap- 
pendix B for details) 

N 

2— ^ J~J g«E<ij> -fij(¥'Ai(n)¥'i3j(")+VBi(n)VA3(n)) 
{<p}=±l n=l 

Xe^i "'(VAi(n)'PAiin+l)+ipBi{n)ipBi{n+l)) ^-[^2) 

where (pM {M = A, B) is the slave-spin in A or B sub- 
lattice since the honeycomb lattice is a bipartite system 
and a = -iln(ef ) > 0. 

Then expanding the above action and working in con- 
tinuum limit, one finds an effective theory as 




where c, r, u are effective parameters depending on mi- 
croscopic details and 1^9^(2^, ''■) ^ fB{x, r) ^ '~p{x, r) since 
only the modes near zero momentum are involved in low- 
energy limit. Therefore, the quantum Ising model in the 
honeycomb lattice can also be described in terms of a t/?^ 
theory in spite of its bipartite feature, which may simplify 
corresponding calculations dramatically. 



B. Hamiltonian of slave- fermion in the honeycomb 
lattice 

Next, we treat the mean-field Hamiltonian Hf. It 
is noted that Hf is a free Hamiltonian when dropping 
the contribution from the constraint. Then the resulting 
Hamiltonian describes a Dirac semimetal which has the 
following formalism in the low-energy limit^l 
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Sf ^ J d^J rf'a:^V^,7^a^V., (14) 

where 70 = h^cTz,!! = cr2 0o-y,72 = h^cr^ with 
I2 the 2x2 identity matrix, a^jCry, being the stan- 
dard Pauh matrix. The Dirac fermion V'a- is defined as 
V'<T = [fi^Jha^ /I<t]^ with the transpose T and we 
also have tpa = V'Jto- /m<t — is the fermion 

\a. M = A^B sublattices near two nonequivalent Dirac 
points located at = -K' = (1, l/\/3)(27r/A/3). It is 
clear that this free Dirac semimetal fixed point is sta- 
ble for any short-range interactions, provided the cou- 
pling parameter is not large. -^ (This stability of free Dirac 
semimetal fixed point can be easily seen from a simple 
RG argument and it does not rely on the existence of 
long ranged Coulomb interactions in Ref. ^2^.) Thus, 
the free Dirac semimetal fixed point with disconnected 
Dirac points can serve as a good starting point in the 
honeycomb lattice at half-filling just like the usual Fermi 
liquid with large Fermi surface. In addition, we should 
emphasize it is the slave-fermion /o- that forms a Dirac 
semimetal but not the physical electrons since in the 
slave-spin approach, the physical electron is a composite 
particle of slave-spin and slave-fermion (cc = fa)- 

Because usually it is more interesting to study the in- 
stability of Dirac semimetal of physical electrons to other 
states via quantum phase transitions, we here assume / 
fermions form a Dirac semimetal whatever the phases of 
slave spins are, and have a sharply defined Dirac quasi- 
particle in the remaining parts of the present paper. (A 
single-site mean-field solution of Eqs. [TU] and [TT] with- 
out the Lagrange multiplier can be found in Appendix D. 
There we find that the f-fermions form a Dirac semimetal 
when slave-spin orders and those fermions decouple from 
each other when slave-spin is disordered. This decou- 
pling is obviously the artifact of the single-site approxi- 
mation. If many-site effect is introduced, we expect that 
the f-fermions form the Dirac semimetal with renormal- 
ized parameters as what has been done in the paper of 
Nandkishore, Metlitski and Senthili^ Even though the 
f-fermions do not form the Dirac semimetal when slave- 
spin is disorder, the orthogonal Dirac semimetal may still 
be considered as a useful intermediate temperature state 
before other long-ranged ordered phases start to form at 
low temperature.) 

C. Nature of physical electrons in the honeycomb 
lattice and orthogonal Dirac semimetal 

To gain some features of the physical electrons in the 
Z2 slave-spin representation for the honeycomb lattice, 
it is helpful to inspect the behavior of the quasiparticle, 
particularly, its single-particle Green's function or equiv- 
alently its spectral function. 

In the case of {t'^) ^ (ordered state of the corre- 
sponding quantum Ising model), it is clear that Cia ~ 



{T^)fiij, which means the slave-fermion corresponds to 
the physical electron when the slave-spin condensates. 
Therefore, the physical electron excitation is a Dirac 
quasiparticle since slave-fermion forms Dirac semimetal 

G{k) = Z'-^ (15) 

where the quasiparticle spectral weight is defined as 
Z = (r^)^ and k'^ — k'^ + uP' with w being the imaginary 
frequency. Obviously, the obtained state with condensed 
slave-spins is just the usual Dirac semimetal since no frac- 
tionalized excitation will appear in physical observable. 
In the language of the gauge theory, the hidden Zi gauge 
field in the Z2 slave-spin representation is confined by the 
Higgs mechanism when the slave-spin condensates, thus 
only gauge singlet of fractionalized particles (/ or r^) are 
allowed in physical excitations due to the confined poten- 
tial generated by the Zi gauge fieldi^i^ One can refer 
to Appendix C for details. 

In contrast, for a vanished (r^) (disordered state of 
the slave spin) and in the low energy limit, the quantum 
Ising model is described by an effective Lp^ theory given 
by Eq. p^ and we obtain the Green's function of slave 
spin as follows^ 

A2+,42+^2 - (16) 

The corresponding spectral function can also be easily 
derived with the form 

A^(k,uj^i^-^)=5{u? ~{^^ ^k^)). (17) 

Clearly, the spectral function of the slave spin has an ex- 
citation gap A. Therefore, the physical electron will also 
acquire a gap with Z — Q. According to the definition of 
the orthogonal metal in the paper of Nandkishore, Metl- 
itski and Senthil)^ if a state has a gap for single-particle 
excitation and the same thermodynamics and transport 
properties as Landau Fermi liquid, it could be identi- 
fied as an orthogonal metal. In our case, the physical c 
electron has an excitation gap while the / fermions form 
Dirac semimetal. More importantly, the / fermions carry 
both charge and spin degrees of freedom of the physical 
c electrons, thus / fermions will contribute to the ther- 
modynamics, charge and spin transports exactly in the 
same way as real electrons. (The contribution of slave 
spins can be neglected in the low energy limit, since they 
are gapped in the disordered state.) Therefore, the c 
electrons are in a gapped semimetal which is named as 
orthogonal Dirac semimetal. 

Evidently, the orthogonal Dirac semimetal will have 
the same thermodynamics and transport properties as 
the usual Dirac semimetal, but it is noted the orthogo- 
nal Dirac semimetal is indeed a Z2 fractionalized state 
with the Z2 gauge field gapped, thus the only active ac- 
tor is the slave-fermion which can be defined as a real 
fractionalized excitation while the usual Dirac semimetal 
is a confined state of Z2 gauge field (see an argument of 
this point in Appendix C). 



6 



Another interestmg point is that, in contrast to the 
usual Dirac semimetal, the Zi fractionalized orthogonal 
Dirac semimetal should have the subtle order, namely, 
the topological orde r^'^i^'* as argued in Ref. [5g. What 
is needed to point out is however that there is a differ- 
ence between our case and the one discussed in Ref. ^56,] 
since the Zi fractionalized orthogonal Dirac semimetal 
has gapless slave-fermion excitations while in Ref. [56| . 
all low-lying excitations for the fractionalized quantum 
spin Hall insulator are gapped. But this distinction seems 
irrelevant to the issue of the topological order and we note 
our case is rather similar to the gapless Z2 spin liquid^^ in 
Ref. Q and the nodal liquid in Ref. [S^l . It is also noted 
that the authors of Ref. [13] argued the effect of matter 
fields to the pure Zi gauge theory and found the topo- 
logical degeneracy still survives even the matter field has 
linearly gapless spectrum (the f-fermions are in just this 
case). The existence of the nontrivial topological order 
means when considering the system defined in a cylin- 
der (torus) , one may find a twofold (fourfold) topological 
degeneracy of the ground state and a topological term 
Stopo — — In 2 in the entanglement entropy J^iiii^ We 
expect such nontrivial topological entanglement entropy 
could be considered as a strong signal to imply the exis- 
tence of orthogonal Dirac semimetal in future numerical 
simulations if one finds a gapped metallic state with no 
symmetry breaking. However, we should emphasize that 
unlike most of slave-particle approaches, here the frac- 
tionalization is irrelevant to the spin-charge separation 
and it seems this feature is general for any fractionalized 
states obtained in Zi slave-spin approach. 



IV. QUANTUM PHASE TRANSITION OF 
ORTHOGONAL DIRAC SEMIMETAL 

Having analyzed the properties of the mean-field the- 
ory, in this section, we proceed to discuss the stability of 
the mean-field treatment, which is crucial for the above 
study of critical behaviors and scaling properties since ef- 
fect of fiuctuation may be fundamentally important near 
QCP. In certain sense, even one finds mean-field decou- 
pling is unstable, it does not imply the disappearance of 
orthogonal Dirac semimetal but only shows the critical 
behavior and phase transition cannot be described by the 
simple theory. 



A. Stability of the mean-field decoupling 

In this subsection, we will study the stability of the 
mean-field treatment which decouples the interacting 
Hamiltonian Eq. into two independent Hamiltonians 
(see, Egs.pU)) and (ITT]) ') with extra self-consistent equa- 
tions. Since the stability of mean-field theory in the Zi 
slave-spin formalism has been discussed in Ref. (36j with 
usual Fermi liquid being the normal state (slave-spin con- 
densates), we will follow their treatment but in our case. 



a crucial difference appears as a Dirac semimetal but not 
the usual Fermi liquid serves as the normal state in the 
honeycomb lattice at half-filling. As a matter of fact, this 
difference leads to an interesting feature that the vanish- 
ing density of states of the Dirac semimetal suppresses 
the effect of Landau damping even in the conditions with- 
out long-range Coulomb interaction, thus the mean-field 
decoupling is more stable than the one in Ref. fSG^ and 
as a result, the critical Fermi surface (Dirac point) is 
obtained rather naturally in the honeycomb lattice, in 
contrast to the case of Ref. [s^ . 

We now show how one can obtain above statement by 
following Ref. [s^. Obviously, the fiuctuation above the 
mean-field theory results from the coupling term between 
the slave-spin and the slave-fermion in Eq. ([3]) 

n^nt = E + '^•cO (18) 

or equivalently one may also use its path integral formal- 
ism 




Since we are only interested in the physics in low-energy 
limit, it is helpful to derive an effective theory as 

5 = y" dTd'^x\L^ +Lf + Lm] (20) 

= (a^^)2 +C2(VV3)2 + V +u^4 (21) 

Lf=Yl V^-TmSmV-. (22) 
Unt = A ^(^2^,^^^^ ^ )^ (23) 

where we have defined a 4 x 4 matrix = 72 -t- 17073 
and L^, Lf come from the free action in Eqs. (jl3p and 
(ITTl) . respectively. And A denotes an effective coupling 
parameter and can be estimated as A ^ too with ao the 
lattice constant. 

It should be emphasized that when the mean-field 
treatment is applicable (In the mean-field theory, one just 
neglects the interacting term Lint-), the quantum criti- 
cal behaviors can only result from the slave spins since 
the corresponding quantum Ising model has a definite 
quantum phase transition while the free Dirac fermions 
contribute no singularity. Therefore, it is more interest- 
ing to study the effect of fluctuations (which comes from 
the interacting term Lint Eq.|23]) to the slave spins and 
a one-loop calculation gives the following correction term 
to action of the slave spins as 

X'J d^+^gf V?(l + ^^)|0(<z)P, (24) 

where 0{x) = ip{x)'^, q = {q,u!) and N = 2 for spin 
degeneracy with c = 1. Then, following Ref. [36|, a 
renormalization-group (RG) argument can be applied in 
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the above action which states the scahng dimension of 
is dini[A2] = | - 1 - with D ^ d + 1, ly being 
the critical exponent in the correlation length ~ Iff I"" 
with g ^ [U — Uc)/t and Uc denoting the critical strength 
of the onsite Coulomb energy at which a QCP locates) of 
the classical Ising models in D. (The critical exponents 
of the quantum Ising model in the space dimension of 
d is identical to the classical Ising model living in d+1 
spatial dimensions and this implies the dynamical critical 
exponent (z) is equal to one.) 

According to the scaling dimension of A^, the correc- 
tion term given by Eq. (1^^ to slave spins, which results 
from integrating over Dirac fermions to lowest order, is 
irrelevant when v > j^p^ is satisfied. For our case, we are 
considering a two-dimensional honeycomb lattice thus 
our system is living in d=2 and one should use v — 0.63, 
a critical exponent for classical three-dimensional Ising 
model. Therefore, it is clear that the induced correction 
term for slave spins is irrelevant in the sense of RG, which 
means the coupling term between the slave-spin and the 
slave-fermion cannot change the low-energy physics ob- 
tained from the simple mean-field approximation. More- 
over, the stability of the mean-field treatment suggests 
that quantum critical behaviors are controlled only by 
the quantum Ising model of slave-spins and all critical 
exponents fall in the Ising universal class. 

However, our result is quite different from the case of 
orthogonal metal^ where generically the mean-field de- 
coupling is only valid in the presence of the long-range 
Coulomb interaction. In our system no such long-range 
interaction is needed. The main reason is that if the 
slave-fermions / form a usual Fermi liquid, they will al- 
ways contribute a standard Landau damping term with 
an extra constant being the density of state (DOS) at 
the Fermi surface for slave-spins. This nonzero DOS 
term is not irrelevant unless one adds the long-range 
Coulomb interaction. In contrast, if the Dirac semimetal 
is considered, as in our case, slave spins will not gain a 
nonzero DOS term but a term given by Eq. ([M)) . just 
like the long-range Coulomb interaction in two space di- 
mensions since for the free Dirac fermions, its DOS at 
the Fermi surface (at the Dirac points in our case) van- 
ishes (In d—2, the Coulomb interaction is V{q) ^ with 

Vix) ~ p^.^^). Moreover, since no Fermi surface (only 
separated Dirac points) appears in our case, we expect 
that no usual Landau damping appears from the long- 
wavelength deformations of the 'Fermi surface' due to 
disappearing of continuous rotation symmetry in lattice 
mo dels <^ 

After all, the mean-field treatment in the last section 
is stable to fiuctuation resulting from coupling between 
the slave-spin and the slave-fermion even without intro- 
ducing the long-range Coulomb interaction. Then we will 
proceed to use the mean-field approximation to analyze 
the critical behaviors and corresponding scaling proper- 
ties near QCP. 



B. Critical behaviors and scaling properties near 
quantum phase transition of orthogonal Dirac 
semimetal 



In this subsection, we will study the critical behaviors 
near the quantum phase transition from orthogonal Dirac 
semimetal to usual Dirac semimetal by using the effective 
action Eqs. (pij) and which refiects the decoupling 
in the mean- field treatment. 

Obviously, the usual Dirac semimetal is well- 
understood and its low-energy excitation is described by 
the Dirac quasiparticle according to the discussion in the 
Sec. mil when approaching the QCP 



G(fc) = Z^, (25) 



where the quasiparticle spectral weight Z — (r^)^ ~ g^^ . 
(One can approach QCP from the orthogonal Dirac 
semimetal with vanished excitation gap as well.) Besides, 
based on a scaling argument in Ref. [68|, the quasipar- 
ticle spectral weight for z = \,d = 2 can be written as 
Z g"^" , where an anomalous dimension rjc for physical 
electrons is defined. Thus, a crucial result in the present 
paper is obtain as rjc — I + rj. (In fact, for d = 2, z = 1, 
this result still holds even if the decoupling approach is 
invalid. To derive above result, we have used two scaling 
relations a + 2f3 + j = 2,v{2 — rj) ~ j and the hyper- 
scaling law 2 — a = v{d + z) with a, /3, 7, rj being the 
critical exponents of the quantum Ising model.) This 
means that at QCP the Dirac quasiparticle is destroyed 
completely {Z = 0) since a large anomalous dimension 
appears {r]c — 1.036,77 ~ 0.036) and the local DOS will 
behave as N{uj) ~ |a;|^"'"'''= which is rather different from 
its counterpart in Dirac semimetal {N{uj) ^ |aj|). In the 
sense of the critical Fermi surface of Senthil,^^ the six 
Dirac points are indeed critical at QCP with nontrivial 
power-law behaviors in the spectral function (or DOS) of 
physical electrons. Thus, we may call them critical Dirac 
points following the similar spirit of the critical Fermi 
surface. However, to our surprised, the quasiparticle pic- 
ture does not break down at QCP because slave spins and 
slave fermions can be considered as the real quasiparticle 
at criticality, respectively. Therefore, quantum critical 
behaviors are readily to be obtained in terms of these 
quasiparticles. For example, at QCP the specific heat is 
C„ ~ which contributes from both slave spins and 
slave fermions. And the electrical conductivity at QCP 
is easily obtained as cr(w) = f f;f which is identical to 
the one in the usual Dirac semimetal since slave fermions 
carry all quantum number of physical electrons and they 
form Dirac semimetal like physical quasiparticles in the 
Dirac semimetal phase. 
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V. EXTENSIONS AND POSSIBLE RELATIONS 
TO OTHER MODELS AND APPROACHES 

In this section, we will briefly argue the instability of 
the orthogonal Dirac semimetal to other phases and dis- 
cussion possible relations of our results obtained in pre- 
vious sections to other models and approaches. A careful 
reader may wonder whether one can construct some ex- 
actly soluble models for the orthogonal Dirac semimetal. 
In fact, it is straightforward to construct possible exact 
soluble models by following the idea of Ref. [s^. The 
only difference needed to be careful is to put their model 
from square lattice to the honeycomb lattice. Here we fo- 
cus on the special features of the honeycomb lattice and 
do not discuss it in detail. 



A. Instability of the orthogonal Dirac semimetal to 
sublattice pairing states and fractionalized quantum 
spin Hall insulator 

In principle, the orthogonal Dirac semimetal proposed 
may be seen as a "Fermi liquid" which can suffer from 
further transitions to other fractionalized states. For ex- 
ample, if a next-nearest-neighbor hopping term is added 
into the Hubbard model, the particle-hole symmetry will 
be violated, thus the terms involving the Lagrange mul- 
tiplier Xi cannot be ignored, which may lead to a paring 
instability of slave fermions fa ■ Particularly, when a pair- 
ing is realized between the next-nearest-neighbor sites, 
the resulting state is like the so-called sublattice pair- 
ing states (SPS)r^'22ii^ which has a gap both for charge 
and spin excitations and has been argued to be a promis- 
ing candidate for gapped spin liquid state found in the 
numerical simulation. However, detail treatment of this 
sublattice pairing states is beyond our present paper and 
it is interesting to attack this issue in future work. 

In addition, one may consider a spin-orbit coupling 
between next-nearest-neighbor sites 

Hso = iX^{{ij))^aa'l^tj<Jaa'4a(^3<T' , (26) 

where A is the spin-orbit coupling and Uij is defined as in 
Ref. [62] with ((•)) denoting next-nearest-neighbor sites. 
Therefore, the resulting model will be the celebrated 
Kane-Mele-Hubbard model which is a good starting point 
for studying interaction effect on two-dimensional topo- 
logical insulators.®^ Obviously, due to this spin-orbit cou- 
pling, the usual Dirac semimetal will smoothly evolve 
into the two-dimensional topological insulator (or quan- 
tum spin Hall state) while the exotic orthogonal Dirac 
semimetal will inevitably open an excitation gap at Dirac 
points, which will give rise to the exotic fractionalized 
quantum spin Hall insulator found by Riiegg and Fietei^ 
This phase could be dubbed an "orthogonal spin Hall 
state" or an "orthogonal helical metal" according to one's 
taste.— Before ending this subsection, we should empha- 
size that a gapped spin liquid has also been found in 
Kane-Mele-Hubbard model when spin-orbit coupling A 



is small ii2r2^ However, our current approach cannot cap- 
ture this subtle issue and we expect future work may 
clarify the nature of this gapped spin liquid and its pos- 
sible relation to the orthogonal spin Hall state. 



B. Comparison the Z2 fractionalization approach 
to the Gross-Neveu theory 

Beside the above possible transitions, it is also noted 
that our results are rather different from the treatment 
of Herbut who did not use fractionalization but started 
with usual order parameters, e.g. spin/charge density 
wave (SDW/CDW), and then coupled them minimally 
to Dirac fermions The resulting theory is described 
by the Gross-Neveu theory with a large anomalous di- 
mension for the order parameter and is argued to be out- 
side the Ginzburg-Landau- Wilson paradigm. But, in the 
effective Gross-Neveu theory, the physical electrons only 
have a small anomalous dimension while our treatment 
gives a rather large anomalous dimension for physical 
electrons since they are fractionalized at QCP and a large 
anomalous dimension is naturally expected. Therefore, 
the large anomalous dimension of physical electrons and 
the related pow-law behavior in the local DOS may be 
considered as a smoke-gun prediction of our Z2 slave-spin 
theory for Hubbard model in the honeycomb lattice when 
comparing to the usual effective Gross-Neveu theory. 



C. Relations to the fermionic dual approach of 
Mott transition and the gauge/gravity duality 

Moreover, it is also interesting to see whether orthog- 
onal Dirac semimetal and the corresponding quantum 
phase transition can be described by extending the dual 
approach developed in Mott transition from Fermi liquid 
to quantum spin liquid.— However, we have tried and 
found that such a description is not feasible due to lack of 
Fermi surface when only half- filling is involved. In con- 
trast, a dual description identical to Refs. [tsI . can 
be constructed when deviating from half-filling since in 
this case the chemical potential is nonzero and a Fermi 
surface instead of Dirac points is formed but this is not 
interesting for our present paper and we will not discuss 
it further. 

The last relation we considered is the link to 
the gauge/gravity duality2Sr22 from string theories. 
Based on some arguments in holographic entanglement 
entropy ) ^^1^^ which is an remarkable application of the 
gauge/gravity duality to calculate entanglement entropy 
of certain dual strongly correlated field theories, only 
gapless metallic states with fermionic dynamical critical 
exponent large than one can be realized in the current 
models in the gauge/gravity duality for no n- Fermi liq- 
uids if a Fermi surface is assumed. Thus, it is clear that 
the exotic orthogonal metal (the simplest non-Fermi liq- 
uid) cannot be found in these models since it is gapped 
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and its fermionic dynamical critical exponent is equal to 
one. For the case of orthogonal Dirac semimetal, to our 
knowledge, no such state has been discovered in exist- 
ing gauge/gravity duality. Therefore, we suspect that 
the nature of exotic orthogonal metal, orthogonal Dirac 
semimetal and their corresponding quantum criticality 
may not be captured by existing models in gauge/gravity 
duality and more study is desired towards this direction. 



VI. DISCUSSION AND CONCLUSION 

In the present paper we have shown that a Zi fraction- 
alized metallic state called orthogonal Dirac semimetal 
can exist in the Zi slave-spin representation of Hubbard 
model in the honeycomb lattice at half-filling. The or- 
thogonal Dirac semimetal can survive when slave spins 
become disordered. This state has the same thermo- 
dynamics and transport as usual Dirac semimetal but 
with gapped singe-partice excitation and nontrivial topo- 
logical order. Moreover, the quantum phase transition 
(QPT) from Dirac semimetal to fractionalized orthogonal 
Dirac semimetal is analyzed by the mean-field decoupling 
and its criticality can be described by the universality 
class of 2-l-lD Ising model. The result that the physi- 
cal electron gains a large anomalous dimension at QCP 
presents the fingerprint of our slave-spin theory. It should 
be emphasized that in contrast to most of slave-particle 
approaches, the Z2 fractionalization involved in the Zi 
slave-spin representation is irrelevant to the spin-charge 
separation. Thus, the orthogonal Dirac semimetal should 
be considered as a simple non-Fermi liquid but without 
spin-charge separation. It seems this feature is general 
for any fractionalized states obtained in Z2 slave-spin ap- 
proach. 

In addition, we have also constructed a path integral 
formalism for the Zi slave-spin representation of Hub- 
bard model and discussed possible relations to effective 
Gross-Neveu theory and gauge/gravity duality. Further- 
more, the instability of orthogonal Dirac semimetal and 
the possible relation to the SPS spin liquid and the frac- 
tionalized quantum spin Hall insulator are briefly ana- 
lyzed as well. In our opinion, the Z^ slave-spin represen- 
tation and its mean-field results seem to work well when 
the onsite V is not too large, otherwise, slave boson or 
slave rotor approach may be more economic to capture 
physics of quantum spin liquids. Besides, we expect the 
exotic orthogonal Dirac semimetal may be realized in the 
sophisticated experiments of ultracold atoms in the hon- 
eycomb optical lattices. ^''"''^ It will be also interesting to 
see whether topological insulator in 3D may provides a 
realization of this state if strong coupling is accomplished 
carefullyi^i^ Moreover, there may also exist an interest- 
ing extension of the present analysis of orthogonal Dirac 
semimetal to the bilayer graphene with a non-vanishing 
density of state.— We hope our findings may be helpful 
for future studies in Zi slave-spin theory and non-Fcrmi 
liquid phases in honeycomb lattice. 
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Appendix A: The explanation of why the slave-spin 
does not carry charge of the physical electron 

Here we should emphasize that although the physical c 
electron has been fractionalized into an auxiliary fermion 
/ct and a slave spin , the quantum number of the elec- 
tron (the spin-i and the charge e) are both carried by 
the / fermion, which is quite different from slave boson 
or slave rotor approaches where the charge and spin de- 
gree of freedom are solely carried by bosonic particles 
and fermionic spinous, respectively. This point is not no- 
ticed until the recent interesting paper— but has crucial 
influence on the correct interpretation of the disordered 
state of the slave spin. As argued in Refs. [36l| , a U(l) 
rotation of physical electron d can only be matched by 
a U(l) rotation of / fermion while the slave spin do 
not change because it is purely real. Therefore electric 
charge must be only carried by / fermion but not the 
slave spin since it corresponds to the Noether charge of 
the U(l) symmetry. 

Appendix B: path integral for the quantum Ising 
model in transverse field 

The quantum Ising model in transverse field is defied 
asi 

ill = -jY.^TlT^3+h.C:)-KY^Tt (Bl) 

where a ferromagnetic coupling J > is assumed and K 
represents the the transverse external field. 

At first glance, one may directly use the coherent state 
of spin operators in constructing the path integral rep- 
resentation, (One can find a brief but useful introduc- 
tion to this issue in Ref. [l|) however, this will lead to 
an extra topological Berry phase term and is not easy 
to utilize practically. An alterative approach is to use 
the eigenstates of spin operator or as the basis for 
calculation;^ One will see this approach is free of the 
topological Berry phase term and give rise to a rather 
simple formalism. Therefore, to construct a useful path 
integral representation, we will follow Ref. (ssf . 

First of all, we consider the orthor-normal basis of Ng- 
Ising spins as 

|a) = |f7i)|f72)|fT2)---|aAr) (B2) 
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with CT, = ±1 and define 



o-j cr) 



T: = |cri;|(T2;|CT3; ■■■\-<j,) ■■■ \(jn) 



(B3) 



(B4) 



Then the partition function Z — Tr{e can be rep- 
resented as 

N 

{t} = ±1 n=l 

where eN=/?. The calculation of (cr(n+l)|e'-^'^' \a{n)) 
is straightforward by exponentiating the rf matrix and 
one gets 

(a{n + l)|e^^^. ^'Win))^ ^{e^^ + e-^VNa,(n + 1)), 



gaCTi(ri)(Ti(n+l)+f) 



(B5) 



where a — ^[lncosh(ei4r) — In sinh(eiir)] and b = 
i[lncosh(eiir) + lnsinh(ei4r)]. Therefore, the resulting 
path integral formalism for the quantum Ising model in 
transverse field is 

N 

1^ ^ "Q geJE(ij> '^.(")'^.(") + Eia<Ti(«)cTi(n+l)+Af,6_ 
{<t} = ±1 n=l 

(B6) 

Further, if one assumes the model is defined in a hyper- 
cubic lattice in space dimension of d, an effective theory 
can be derived as 



where ^ = [^Y^"" with oq being the lattice constant 



and 



2 _ Joq 



Moreover, in the effective theory. 



corresponds to while gives the kinetic energy term 
in imaginary time. Then, the standard theory is ob- 
tained by relaxing the hard constraint 0^ = 1 while in- 
troducing a potential energy term. 



Z = \ D4>e 



- /(ird''x[(a^0)2+c^(V0)^+r0^+u0*] 



, (B8) 



where r, u are effective parameters depending on micro- 
scopic details. 



Appendix C: Z2 gauge theory formalism for Z2 
spin-slave representation 



the slave-spin and the slave-fermion in Eq. one ob- 
tains 

J dTY^[f,„^rf^a+a{^T<P^)^] 

i 

- / dr'^iJijlfiipj +iijficrf]cr +C.C.) 

+ j dT-Y..hfU, (CI) 

where the summation in the imaginary time is trans- 
formed to an integral and J^j, tij are the auxiliary field 
introduced in the Hubbard-Stratonovich transformation, 
respectively. The usual mean-field approximation can be 
recovered by treating Jij,iij as static variable. Beyond 
mean-field approximation, one may allow a phase fluctu- 
ation in the static , Uj , which are labeled as J.^f^ , tfj^ 
to emphasize they are solutions of the mean-field approx- 
imation. However, since the original Hamiltonian Eq. ^ 
only has the local Z2 symmetry, the phase fluctuation is 
nothing but the expected Z2 gauge fleld deflncd in the 
link between two siteSf^ Therefore, we have 



IMF 



(C2) 



and 



where aij denotes the dynamical Z2 gauge field and terms 
without dynamics are neglected. Much physics can be 
seen from the above action. If the slave-spin {ip) conden- 
sates, the Z2 gauge fleld will be conflned by the Higgs 
mechanism while the slave-spin and slave-fermion have 
to bind into physical electrons as the only low-energy 
excitation. In contrast, the disordered slave-spin ((^) is 
gapped and can be integrated out to generate an effective 
"kinetic energy" term for the Z2 gauge fleldi^ In prin- 
ciple, the Z2 gauge fleld can be in the deconflned state 
in 2-1- ID if the effective kinetic energy is large enough to 
suppress the fluctuation of the Z2 gauge fleld. Then, in 
the deconflned state, one can treat the gauge fleld as a 
static background (The excitation of the gauge fleld, the 
Z2-vortex (vison) is gapped.), where the slave-fermion is 
almost free and can be deflned as real excitation, thus 
one has a Z2 fractionalized state indeed. 



Here, we derive the Z2 gauge theory formalism for Z2 
spin-slave representation by using Eqs. ([5]) and ©• 

First, it is noted in Eq. ^ the constraint can be 
dropped because we are only interested in non-magnetic 
solutions. Then using the familiar Hubbard-Stratonovich 
transformation to decouple the coupling term between 



Appendix D: A self-consistent mean- field solution 

In this Appendix, we give the single-site self-consistent 
mean-fleld solution for Eqs. [lU] and [11] without the 
Lagrange multiplier. 
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Following the treatment of Sec. 4 in Ref.fs^ and using 
Eqs. [To] and [11], the self-consistent equation for g = 
(r^)^ is obtained as 

4e 

1 = = (Dl) 

V(?//4)2 + We^g 

where we have defined e = dep{e)e with p{e) denoting 
density of states of free f-fermions. When g = 0, one finds 



the critical coupling Uc = 16e ~ 12.444t. For U > Uc, 
the slave-spin does not form ferromagnetic order while 
the ferromagnetic order appears with nonzero g = 1 — 
(U/Uc)'^ ioT U < Uc- As for the f-fermions, they can be 
described by the following Hamiltonian 

H = -tgJ2{flf,.+H.c.). (D2) 
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